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ABSTRACT 


Context. The 16 Cygni system is composed of two solar analogues with similar masses and ages. A red dwarf is in orbit around 16 
Cygni A, and 16 Cygni B hosts a giant planet. The abundances of heavy elements are similar in the two stars, but lithium is much 
more depleted in 16 Cygni B than in 16 Cygni A, by a factor of at least 4.7. 

Aims. The interest of studying the 16 Cygni system is that the two star have the same age and the same initial composition. The 
differences currently observed must be due to their different evolution, related to the fact that one of them hosts a planet while the 
other does not. 

Methods. We computed models of the two stars that precisely fit the observed seismic frequencies. We used the Toulouse Geneva 
Evolution Code (TGEC), which includes complete atomic diffusion (including radiative accelerations). We compared the predicted 
surface abundances with the spectroscopic observations and confirm that another mixing process is needed. We then included the 
effect of accretion-induced fingering convection. 

Results. The accretion of planetary matter does not change the metal abundances but leads to lithium destruction, which depends 
upon the accreted mass. A fraction of the Earth’s mass is enough to explain the lithium surface abundances of 16 Cygni B. We also 
checked the beryllium abundances. 

Conclusions. In the case of accretion of heavy matter onto stellar surfaces, the accreted heavy elements do not remain in the outer 
convective zones, but are mixed downwards by fingering convection induced by the unstable p-gradient. Depending on the accreted 
mass, this mixing process may transport lithium down to its nuclear destruction layers and lead to an extra lithium depletion at the 
surface. A fraction of the Earth's mass is enough to explain a lithium ratio of 4.7 in the 16 Cygni system. In this case beryllium is not 
destroyed. Such a process may be frequent in planet-hosting stars and should be studied in other cases in the future. 


1. Introduction 


The bright solar analogues 16 Cygni A (HD 186408, HR 7503) 
and 16 Cygni B (HD 186427, HR 7504) for many reasons rep¬ 
resent a very interesting stellar system. While a red dwarf, 16 
Cygni C, is in orbit around the first c omponent 16 Cygni A 
iTurner et alj200lilPatience et alJ2002h . the second component, 
16 Cygni B, hosts a giant Jovian planet with minimum mass of 
1.5 Mj up located on an eccentric_orbit_(e=0.63), with an orbital 
period of 800.8 days dCochran et al.ll 1997h . The two main stars 
are separated widely enough, with an orbital period longer than 
18,000 years (lHauser & Marcvl[l999t) . to be studied in the same 
way as two isolated stars, with no common dynamical effects. 

This situation allows for precise differential studies between 
a planet-hosting star and a non-planet-hosting star with similar 
birth conditions. The red dwarf around 16 Cygni A may be the 
reason why no accretion disk has developed around it, whereas 
a planetary disk remained around 16 Cygni B, including the ob¬ 
served giant planet, and probably smaller as yet unobserved bod- 


The two main stars of the 16 Cygni system have been stud¬ 
ied in many ways, using spectroscopy, interferometry, and as- 
teroseismology. The abundances of the heavy elements in these 
two stars are ve ry similar. Although a v ery small difference has 
been claimed bvlRamifez et al.ll2.01 ll and lTucci Maia et akll20l4 
ISchuler et alJ (1201 ll) found them to be indistinguishable. On the 


other hand, the surface lithium abundance of 16 Cygni B is lower 
than that of 16 Cygni A by at least a factor 4.7 (iKing et al.ll997l) . 
These observations lead to several open questions, which remain 
to be answered for a better understanding of these stars. 

The interest of this study is that these stars have the same 
birth site and the same age, with masses of the same order, so 
that their past evolution is similar for most aspects. The observed 
differences between them must basically be due to the planetary 
disk around B. For this reason, the detailed study of this stellar 
system helps understanding the differences between stars with 
and without planetary disks. 

The present paper is motivated by two considerations. The 
first one is that none of the previous modelling of these two stars 
considered atomic diffusion including the radiative acceleration 
on each element. Most stellar evolution codes currently include 
the atomic diffusion of helium (without radiative accelerations), 
some of them also include it for heavy elements, but very few 
consider the radiative accelerations. 

The second consideration refers to the consequences of the 
accretion of heavy matter onto stars, which may occur when the 
star has a planetary disk. Many studies and past publications 
assumed that the accreted matter remains inside the outer stel¬ 
lar convective zone, so that accretion can lead to an increase of 
the heavy element abundances, as well as to an increase of the 
lithium abundance. This assumption is not valid, as shown i n 
detail by IVauclarH (120041) . iGaraudl (1201 ll) . iTheado et al.l (1 201 2l) . 
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and iDeal et al.l ( 20138) . When heavy matter falls onto the star, it 
creates an inverse gradient of molecular weight, which leads to 
a double-diffusive instability now called fingering (or thermoha¬ 
line) convection. The heavy elements are mixed downwards until 
the mean molecular weight gradient becomes nearly flat. In most 
cases, no signature of the accreted heavy elements re mains at the 
surfac e. Meanwhile, as computed in detail bv lTheado & Vauclairl 
d2012l) , the induced mixing may lead to an extra lithium deple¬ 
tion in the star. As a consequence, the accretion of heavy matter 
cannot lead to any increase of lithium at the surface of a star, but 
may conversely lead to a decrease of its observed abundance. 

The observational results on the two main stars 16 Cygni 
A and B are presented in Sect. 2 together with a discussion of 
the still-unsolved questions. In Sect. 3 we present models of 
16 Cygni A and B that fit the observed seismic frequencies. 
We compare the parameters of these models with previously 
published ones. The surface abundances of heavy elements and 
lithium that are obtained after diffusion in these models are dis¬ 
cussed in Sect. 4 . Finally, in Sect. 5 we show that the accretion 
of metal rich planetary matter at the beginning of the main se¬ 
quence on 16 Cygni B may explain the lithium difference be¬ 
tween the two stars. A summary and discussion of all these re¬ 
sults are given in Sect. 6. 

2. Observational constraints 

The position in the sky of the 16 Cygni system allowed seis¬ 
mic observations with the Kepler satellite. More than 40 modes 
of degree 1=0, 1, 2, and 3 could be detected. Analyses with the 
Asteroseismic Modelling Portal (AMP) and comparisons with 
other seismic studies led to p recise values of the masses, radii, 
and ages of the two stars dMetcalfe et al.ll2012 ). Moreover, the 
seismic observ atio ns allowed measurements of their stellar rota¬ 
tion periods dDavies et all2015h . 

The effective temperatures and gravities of the tw o stars were 
derived from spectroscopic observations (e.g. IRamfrez et aTl 

1201 lb ISchuler et akll201 it pTucci Maia et al.ll2014l) . We must in¬ 
sist, however, that asteroseismology leads to a log g value with 
a much better precision than spectroscopy. As we show in Sect. 

3, all the models that closely fit the observed seismic frequen¬ 
cies, even if they have different masses, radii, helium abun¬ 
dances, ages, and luminosities, have the same log g value with 
a precisi on of 0.01. T he determinations of the bo lometric mag- 
nitudes ( Torresi12010 ) and Hipparcos parallaxes (ivan Leeuwenl 
1 20071) lead to precise values of the luminosities of the two stars. 
They are also bright enough for their_radii to be determined by 
interferometric techniques (White et al.>li2013 ). The results are 
given in Table [I] 

Detailed determinations of their element abundances have 
been given by several authors. From spectroscopy with high 
signal-to-noise ratio of the 10 m K eck 1 telescope and HIRES 
echelle spectrogranh. lSchuler et al.l (1201 lh determined the abun¬ 
dance of 15 heavy elements in both stars and found them in¬ 
distinguishable. Onjhe other hand, iRamfrez et al.l (1201 lh and 
ITucci Maia et al.l (12014l) claimed that 16 Cygni A is slightly 
more metal rich than 16 Cygni B based on spectra with high 
resolution and high signal-to-noise ratio that were obtained with 
the R.G. Hull coude spectrograph on the 2.7m Harlan Smith tele¬ 
scope at Me Donald Observatory. 

A more striking difference between the two stars, confirmed 
by all spectroscopic observations, is that the star hosting a planet, 
16 Cygni B, has an abundanc e of lithium at least four times lowe r 
than the star without a planet dFriel et al.ll993l:lKing et al.l 19971) . 
While both stars are lithium depleted compared to F stars and to 


Table 1. Properties of 16 Cygni A and B from the literature 



16 Cygni A 

16 Cygni B 

Teff(K) 

5825 ±50" 

5750 ± 50" 


5813 + 18* 

5749 ±17* 


5796 ± 34 e 

5753 ± 30" 


5839 ± 42 d 

5809 ± 39 rf 


5830 ± l f 

5751 + 6/ 

logg 

4.33 ± 0.07" 

4.34 ± 0.07" 


4.282 ± 0.017* 

4.328 + 0.017* 


4.38 ± 0.12" 

4.40 + 0.12" 


4.30 ± 0.02-/ 

4.35 ± 0.02 f 

[Fe/H] 

0.096 ± 0.026" 

0.052 + 0.021" 


0.104 + 0.012* 

0.061+0.011* 


0.07 ± 0.05" 

0.05 ± 0.05" 


0.101 ±0.008/ 

0.054 ± 0.008 f 

A(Li) 

1.27 + 0.05' 

< 0.6' 

A(Be) 

0.99 + 0.08/ 

1.06 + 0.08/ 

Mass (M 0 ) 

1.05 + 0.02 s 

1.00 + 0.01* 


1.07 + 0.05'/ 

1.05 + 0.04'/ 


1.11 ±0.02* 

1.07 + 0.02* 

Radius (R 0 ) 

1.218 + 0.012'/ 

1.098 + 0.010'/ 


1.22 + 0.02" 

1.12 + 0.02" 


1.243 + 0.008* 

1.127 + 0.007* 

Luminosity (L e ) 

1.56 ± 0.05 s 

1.27 + 0.04* 

Age (Gyr) 

7 i c+0.04fr 

7 ‘ 1 -1.03 

7 o/:+0.69fr 

7 -0.33 


6.9 ±0.3* 

6.7 ± 0.4* 

z, 

0.024 ± 0.002* 

0.023 ± 0.002* 

Yi 

0.25 + 0.01* 

0.25 + 0.01* 

v sin i (km.s -1 ) 

2.23 ± 0.07* 

1.27 + 0.04* 

P rot (days) 

23.8+J f 

23.2 4 “ 2 5 * 

Planet detected 

no 

yes* 


No tes. (n) IRamfrez et all d2009|) : (b> IRam frez et all d201 lh ; 
<c) ISchuler et all 1201 ll) : <" > IWhite et all (120131) . seismic de- 
termination ; (e) IWhite et all (120131) . in terferometric d etermi na¬ 
tio n : w ITucci M aia et all lioi ^ : iMetc alfe et al.l d2012|) ; 

w iDavies et a 1J d2QI5l) ; (l) iKing et all d!997l) : w iDelivannis et al.l 
d2000i) : m ICochran et all dl997l) 


the meteoritic value, detailed measurements show that 16 Cygni 
A is slightly less depleted than the Sun (log N(Li) - 1.27 com¬ 
pared to 1.05 for the Sun), whereas 16 Cygni B is mor e depleted 
(log N(Li) < 0.60). On the other hand. IDelivannis et al.l (120001) 
found that the Be and B abundances are the same in the two stars 
within the limits of the uncertainties. 


3. Asteroseismic studies and stellar models 
including radiative accelerations 

3.1. Stellar models 

We used the Toulouse Geneva Evolution Code (TGEC) to 
compute stellar models that fit the seismic observations of 16 
Cygni A and B. This code performs complete computations 
of atomic diffusion, including radiative accelerations, for 21 
species, namely 12 elements and their main isotopes: H, 3 He, 
4 He, 6 Li, 7 Li, 9 Be, 10 B, 12 C, 13 C, 14 N, 15 N, le O, 17 0, 18 0, 20 Ne, 
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log(A M/M,) 


Fig. 1 . Profiles of the radiative accelerations of three elements Fe, Mg, 
and C as a function of the mass fraction in 16 Cygni A. 


22 Ne, 24 Mg, 25 Mg, 26 Mg, 40 Ca and 56 Fe ( Theado et aljl20H . 
The diffusion coefficients used in the code are those derived by 
iPaauette et al.l (119861) . 

The Rosseland opacities are recalculated inside the model, 
at each ti me ste p and at every mesh point, using OPCD v3.3 and 
data from ISeatonl ( 120051 ) . to take into account the local chemi¬ 
cal composition. In this way, the stellar structure is consistently 
computed all along the evolutionary tracks, as well as the indi¬ 
vidual radiative accelerations of C, N, O, Ne, Mg, Ca, and Fe. 
This is done by using the improved semi-a nalytical prescrip¬ 
tion proposed by lAlecian & Lebla nc ( 20041 ) ■ A more detailed 
discussion of these computations is given in iDeal et all (l2015h 
(in prep). 

The equation of state used in the code is the OPAL2001 equa¬ 
tion ([Rogers & Navfonovll2002|). The nuclea r reaction rates are 
from the NACRE compilation (lAnguloll999l) . The mixing length 
formalism is used for the convective zones with a mixing length 
parameter of 1.8, as needed to reproduce solar models. 

As an example. Fig. Q] displays the g ra d profiles for C, Mg, 
and Fe compared to gravity for one of the computed models of 16 
Cygni A. This model best fits the asteroseismic observations, as 
we discuss in Sect. 3.2. The dashed vertical line represents the 
position of the surface convective zone. Below the convective 
zone g ra d is at least 1/3 smaller than g, so that the effect of the 
radiative accelerations is small and generally negligible inside 
these stars. 

3.2. Best models from asteroseismic fits 

The stellar oscillation_modes_ were derived for each model us¬ 
ing the PULSE code (iBrassard & Charpinetl2008l) . The frequen¬ 
cies were corrected for surface effects in the way proposed by 
iKieldsen et al.l d2008l) . 

We computed models with masses ranging from 1.05 to 
1.14M 0 . The initial helium mass fraction Y, was varied from 
0.245 to 0.26 and the heavy element mass fraction Z, from 0.023 
to 0.025 (initial mass fraction of all elements heavier than he¬ 
lium). 

Evolutionary tracks computed for two different initial com¬ 
positions Yj- 0.25; Z,=0.024 and T,=0.26; Z,=0.024 are pre¬ 
sented in Fig. |2] 


The observational uncertainties on the oscillation frequen¬ 
cies for l = 0 to l = 2 modes lie between 0.1 to 1.45 pHz 
(jMetcalfe et al.ll2012l) . The derived large separations Av for 16 
Cygni A and for 16 Cygni B are Av„/, V = 103.56 + 0.10 jiHz. 
and Av 0 b s = 117.17 ± 0.10 /iHz, respectively. The models with 
large separations that fit these values are represented by blue 
dots in Fig. [2] They include the corrections for surface effects as 
proposed bv Tkieldsen et al.l d2008l) . Neglecting this effect would 
lead to models slightly above the blue dots, with a small differ¬ 
ence in age of 150 million years on average. We note that all the 
models presented in Fig. [2] lie on horizontal lines, which means 
that they all have about the same surface gravity. This is a well- 
known result of asteros eismology. The asymptotic treatment of 
the oscillations (e.g. lTassoulll 19801) shows that the large separa¬ 
tion Av directly gives the average stellar density. In the range 
of our possible models, the variations in radii are small so that 
log g is also nearly constant. In other words, if any parameter of 
the model is modified, for instance, the initial Y value, the other 
parameters such as age adjust to obtain a model with the same 
large separation, and thus the same gravity. Spectroscopy is only 
used in our description to constrain the effective temperatures. 

The uncertainty on the log g values derived from the ob¬ 
served seismic frequencies for each track is about 10 2 . This is 
the uncertainty on the position of the models represented by blue 
dots in Fig. [2] The corresponding error bars would lie inside the 
printed symbols. We note that th e seis mic log g values lie out¬ 
side the ranges given bv lRamfrez et al.l (1201 ll) (see Fig. [2]), which 
suggests that their uncertainties are underestimated. 

We then derived the best of all these models for both stars. 
We first did it independently for each Y, value. The best models 
(represented by squares) were selected first because their small 
separations bvo ,2 best fit the observed ones, second because their 
echelle diagrams best fit the observed ones, according to^ 2 min¬ 
imisations performed between the observed and the modelled 
frequencies. 

When compared with spectroscopic observations (see Fig. 
0, it is clear that the best m odels for Y,=0.25 lie inside the 
box reported bv ISchuler et al.l (12011 *) for 16 Cygni A alone, not 
for 16 Cygni B, and th at they are outside the box derived by 
iTucci Maia et alJ d2014l) for both stars. In contrast, for T,=0.26 
lie inside all the boxes for both stars. For this reason, we find that 
the 17=0.26 value is more probable than that of 17=0.25. 

The echelle diagrams corresponding to the best case com¬ 
puted with K,=0.26 are presented in Fig. [3] for 16 Cygni A 
(left panel) and 16 Cygni B (right panel). The fits between 
the computed and observed frequencies are very good for both 
stars. The corresponding stellar masses are Ma - 1.10M o and 
M b = 1.06M o . 

These values and the other parameters obtained for these best 
models are given in Table 0 The uncertainties we list in this 
table correspond to the possible range of parameters of all the 
computed models, which have large and small separations in the 
observational uncertainties. 

We must note that the surface helium abundances derived 
from these computations are slightly subsolar, while the spec¬ 
troscopic parameters given in the literature have been computed 
using a solar helium value. It would be interesting in the future to 
iterate with spectroscopists and see how such a helium difference 
could influence the derived effective temperature range. 
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Fig. 2. Evolutionary tracks for models from 1.05 to 1.14 Af 0 (from right to left) with Z, = 0.024 and F, = 0.25 (left panel) and evolutionar y tracks 
for mo dels from 1.05 to 1 .12 M Q (from right to left) with Z, = 0.024 and F, = 0,26 (right pa nel) . The error boxes ar e those of|Ramrre zetaJ 
d2009h (red dashed lines). IRamfrez et akl d20 1 1 1 ) (green dot-dashed lines). [Schuler et ah (1201 ll) (blue dotted lines) and lTucci Maia et alj <2014 1 
(black dotted lines). The blue dots indicate models with the right large separation, taking Kieldsen et al. (2008) corrections into account. The blue 
squares correspond to models that also have the right small separations and best fit the Eche lle diag ram. The black thick segments of each line 
indicate the models whose radii are consistent with the interferometric determinations of lWhite et alj ( 2013 ). 
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Fig. 3. Echelle diagrams for 16 Cygni A (left panel) and 16 Cygni B (right panel). The observed frequencies are represented by crosses.The 
frequencies computed for the models with Y=0.26 and Z=0.024 are represented by blue dots (1=0), green triangles (1=1), red squares (1=2), and 
black diamonds (1=3). 




4. Heavy elements and lithium abundances in 16 
Cygni A and B, observational and theoretical 
discussions 

4.1. Spectroscopic observations 

The abundances of heavy elements in 16 C ygni A and B have 
been a subject of debate for 15 years. iDelivannis et al.l ( 2000 ) 
reported an iron overabundance higher in 16 Cygni B than in 16 
Cygni A, which has_not been confirmed by more recent papers. 
Whereas ISchuler et al.l (1201 ll) did not report any difference be- 
tween 16 Cygni A an d B for the abundance s of he avy elements, 
IRamfrez et al.l (1201 ll) and iTucci Maia et al.l (12014l) claimed that 
the heavy elements are overabundant in 16 Cygni A compared to 
16 Cygni B. In any case, all observers agree that lithium is more 


depleted in 16 Cygni B than in 16 Cygni A by a large factor. 
Lithium has not been detected in B, whereas its abundance in A 
is slightly higher than that of the Sun. This difference between 
these very similar stars is difficult to account for using traditional 
explanations of lithium depletion in G stars. 


In the following, we discuss the influence of atomic diffusion 
and mixing on the lithium and heavy elements abundances, and 
we show how the lithium difference between the two stars may 
be explained. We also discuss the consequences of these effects 
for heavy elements. 
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Table 2. Properties of 16 Cygni A and B from this work 



16 Cygni A 

16 Cygni B 

Teff(K) 

5821 +25 

5747 ± 25 

log,? 

4.29 + 0.01 

4.36 + 0.01 

Mass (Af e ) 

1.10 + 0.01 

1.06 + 0.01 

Radius ( R 0 ) 

1.24 + 0.01 

1.13+0.01 

Luminosity (L Q ) 

1.58 + 0.03 

1.25+0.03 

Age (Gyr) 

6.4 ±0.4 

6.4 ± 0.4 

z, 

0.024 

0.024 

Y, 

0.26 

0.26 

7 a 
^ surf 

0.0221 

0.0223 

V a 

* surf 

0.2226 

0.2265 


Notes. Values at the age of best models 



Fig. 4. Abundance evolution of C, Mg, and Fe in our best models of 
16 Cygni A and B, under the influence of atomic diffusion below the 
convective zone, when no mixing is taken into account. 


4.2. Influence of atomic diffusion on the abundances of 
heavy elements and lithium 



log(AM/Mj 

Fig. 5. Abundance profiles C, Mg, and Fe inside the stars under the 
same conditions as for Fig. [4] 



Fig. 6. Abundance profiles of lithium under the same conditions as for 

Fig-El 


Our models include the computation of detailed atomic diffu¬ 
sion for a large number of elements, as discussed in Sect. 3.1. 
Figure [4] presents the surface abundance evolution for three of 
these elements, C, Mg, and Fe, in both 16 Cygni A and B. As 
the mass of 16 Cygni A is higher than that of 16 Cygni B, the 
convective zone is smaller and the diffusion processes faster. As 
a consequence, the surface abundances predicted by the mod¬ 
els are higher in 16 Cygni B than in 16 Cygni A by ~ 0.002 
dex at the age of our best models (vertical line). Later on in the 
evolution, at around 8 Gyr, the situation is reversed because the 
surface convective zone sinks more quickly in the more massive 
star. This occurs in models that are too old to account for the 
observations, however. 

Figure 0 displays the abundance profiles for the same three 
elements C, Mg, and Fe in the best models for 16 Cygni A and 
B. The bottoms of the convective zones are shown as vertical 
dashed lines. The difference of heavy elements abundances is 
due to the difference in surface convective zone depth, which is 
deeper in 16 Cygni B than A (grey dashed line). Atomic diffu¬ 


sion coupled with nuclear destruction also has an effect on the 
abundance of 7 Li (see Fig. [6]). In this case, the final abundance 
of 7 Li is lower in 16 Cygni B than in 16 Cygni A. This is not 
enough to account for the observations. First lithium is depleted 
in 16 Cygni A by a factor smaller than two, while in the real star 
it is depleted by 100. Second the depletion ratio between the two 
star is only ~ 1.4, while the observations show a ratio of at least 
4.7. Extra mixing below the convective zones is clearly needed 
to account for the observations. 


4.3. Mixing processes below the outer convective zone 

Many papers have been published in the past 30 years on the sub¬ 
ject of the mixing processes that may occur below the convective 
zones of solar type stars and result in lithium destruction in their 
out er layers. The most common one is rotation-induced mixing 
(e.g Vauclairf 1 988l Pinso nneault etall 1 1 99(1 ICharbonnel et all 
1 1 9921 1 1 994ICastro et al Jl2009h . 
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The rotation periods of 16 Cygni A and B have recently 
been measur ed _with high precision using asteroseismology 
dDavies et al.ll2015h . They are very similar, 23.8 days for A, 23.2 
days for B, with a difference smaller than the uncertainties. The 
authors have also determined their inclination angles, 56 degrees 
for A, 36 degree for B. Using the interferometric radius, the re¬ 
sulting v sin i is 2.23 km.s~' for A and 1.27 km.s -1 for B. 

When studying rotation-induced mixing, it appears that the 
efficiency of the mixing is related to the local linear rotation ve¬ 
locity in such a way that it increases with radius. The bottom 
of the convective zone lies at a larger radius in A than in B, 
but the mixing is also more efficient. Simple expressions of the 
mixing diffusion coefficients (e.g. lZahnill992l) include the factor 
(Q 2 R 3 /GM), which is about 1.3 larg er for 16 Cygni A than for 
16 Cygni B. As already discussed bv lDelivannis et aD d2000l) . it 
does not seem possible to account for the lithium abundance ra¬ 
tio between A and B by such a simple process, ft would not be 
realistic to invoke a stronger rotation-induced mixing effect in B 
than in A simply to account for the observations. Another pro¬ 
cess is clearly needed. We do not intend to discuss all these pro¬ 
cesses here. We only simulate turbulent mixing by using a turbu¬ 
lent diffusion coefficient adjusted to obtain a lithium destruction 
by a factor 100 in 16 Cygni A, and we used the same to deduce 
the lithium destruction in 16 Cygni B. As in the simulations of 
iRicher etakl (l2000i) , we used a simple form Dj - toD( H e ) (jf) 
with n = 3, a) — 325, Z)(He) the helium diffusion coefficient and 
po the density at the bottom of the outer convective zone. 

Figure [7] displays the resulting surface lithium abundance 
evolution in the two stars (green solid line and blue dashed line). 
As expected, lithium is more destroyed in 16 Cygni B than in 16 
Cygni A, but only by a factor 2.9, not enough to account for the 
observations of an abundance ratio higher than 4.7 at the age of 
the stars. 

Meanwhile, the abundance of heavy elements is also modi¬ 
fied, but the abundances in 16 Cygni B may not become lower 
than in 16 Cygni A (see Fig. ©. 

5. Fingering convection and lithium destruction 
induced by planetary accretion 

We now return to the fact that the star 16 Cygni B hosts a giant 
planet while 16 Cygni A does not. As discussed in the introduc¬ 
tion, the red dwarf that orbits 16 Cygni A may be the reason for 
the fact that the main star could not develop any planetary disk. 

In this situation, it is quite possible that in its early period on 
the main sequence, 16 Cygni B was able to accrete some mat¬ 
ter from its planetar y disk, which could not oc cur for 16 Cygni 
A. As discussed by (iTheado & Vauclairll2012l) . the accretion of 
heavy planetary matter onto a stellar surface builds an unstable 
compositional gradient at the bottom of the surface convective 
zone that triggers fingering (thermohaline) convection. This mix¬ 
ing leads to complete dilution of the heavy matter inside the star, 
so that no signature appears at the surface for the heavy elements. 
On the other hand, it may lead to extra lithium depletion because 
the mixing continues down to the destruction layers. 

5.1. Computations of fingering convection 

Fingering convection may occur in stars every time a local ac¬ 
cumulation of heavy elements appears in the presence of a sta¬ 
ble temperature gradient, ft occu rs in particular in the case of 
the accretion of planetary matter (IVauclairll2004t lGaraudll201 It 
iDeal et al .1120131) . 



Fig. 7. Evolution of lithium surface abundances including a mixing at 
the bottom of the outer convective zone needed to reproduce the 16 
Cygni A abundance in the two models (green solid line and blue dashed 
line) (Sect. 4.3). Evolution of lithium surface abundances including ac¬ 
cretion of 0.66 M s for 16 Cygni B (blue dotted line) at the beginning 
of the main sequence, but with the same mixing as befor e (Sect. 5). 
Abundances determined using observations from lKing et al .’ < 1997 1 are 
represented by black crosses. 



Fig. 8. Abundance evolution of C, Mg, and Fe in our best models of 
16 Cygni A and B, under the influence of atomic diffusion below the 
convective zone, when mixing is taken into account. 


Fingering convection is characterised by the so-called den¬ 
sity ratio Rq , which is the ratio between the thermal and p- 
gradients: 


Ro 


V~V fl d 

V, 


The instability can only develop if this ratio is higher than one 
and lower than the Lewis number, which is the ratio of the ther¬ 
mal to the molecular diffusivities. In this case, a heavy blob of 
fluid falls down inside the star and continues to fall because it ex¬ 
changes heat more quickly than particles with the surroundings. 
If Rq is smaller than one, the region is dynamically convective 
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(Ledoux criterium), and if it is larger than the Lewis number, the 
region is stable. 

Various analytical treatments of fingering convection in stars 
were given in the past, leading to mixing coefficients that dif¬ 
fered by orders of magnitude (lUlrichlIl972t iKippenhahn et al.1 
Il980h . More recently, 2D and 3D numerical simulations 
were performed, converging on coeffi cients of similar orders 
(lDenissenkovll2010i : lTraxler et al.11201 ll) . 

Here we used the recent prescription given by iBrown et al.l 
( 20131). which has be en confirmed by the 3D simulations of 
IZemskova et al.ll20l4 

5.2. Effect of accretion-induced fingering convection on the 
lithium abundance of 16 Cygni B 

We computed stellar models of 16 Cygni B with the assumption 
of accretion of planetary matter at the beginning of the main- 
sequence phase. These models included the treatment of finger¬ 
ing convection in the case of an inversion of the mean molecular 
weight gradient. We tested different accretion masses to char¬ 
acterise their effect on the lithium surface abundance due to 
fingering convection. In these computations, the accreted mat¬ 
ter was assum ed to have an Earth-like chemical composition 
(lAllegre et al.lll995l) . Changing the relative abundances of the 
accreted heavy elements has a very weak effect because their 
contribution to the mean molecular weight may be slightly dif¬ 
ferent. For example, decreasing the iron abundance by a factor 
two compared to other elements modifies the computed accre¬ 
tion mass by less than 10%. 

Abundance profiles of lithium after the accretion of various 
masses of planetary matter are presented in Fig. [9] 

An accreted mass lower than or equal to 0.6 M® does not 
have any real effect on the lithium abundance because it reduces 
its surface value by a factor lower than 10 percent. The mixing is 
not efficient enough to mix the lithium down to the destruction 
layers. For masses higher than 0.6 M®, the lithium destruction 
becomes important and reaches a factor 3.5 for an accreted mass 
of 0.66 M® and more than a factor 100 for an accreted mass of 
1 M®. 

Thus a small amount of planetary matter is enough to sig¬ 
nificantly reduce the lithium surface abundance in this type of 
stars. 

After the accretion episode in 16 Cygni B, the lithium abun¬ 
dance continues to decrease due to other extra mixing processes 
in the two stars, as discussed in the previous section. We plot in 
Fig- E] the lithium surface abundance evolution with the assump¬ 
tion of an accretion episode of 0.66 M® (blue dotted line), which 
is enough to account for the observations. As the lithium obser¬ 
vation in 16 Cygni B is only an upper limit, any accretion mass 
higher than this value can explain the abundance difference in 
the two stars. 


5.3. Beryllium 

Bery llium has been detected in 16 Cygni A and B by 
iDelivannis et al.l ( 2000 ). who found that the difference between 
the two abundances, if any, must be smaller than 0.2 dex. Beryl¬ 
lium is destroyed by nuclear reactions at a temperature of T ~ 
3.5 10 6 /C, which is higher than the temperature needed to de¬ 
stroy lithium. The fact that it is not depleted gives an upper limit 
on the accreted mass that must not lead to mixing down to its 
destruction region. 



Fig. 9. Lithium abundance profiles after the accretion of different 
masses in the model of 16 Cygni B: no accretion (blue solid line), 
0.6 M e (green dashed line), 0.66 M 0 (red dotted line), and 1 M e (cyan 
dotted-dashed line) at the beginning of the MS. 


To characterise this effect, we plot the same figure as Fig. [9] 
for beryllium (see Fig.fTOli. An important result is that for an ac¬ 
cretion of 0.66 M® (red dotted line) beryllium is not destroyed by 
nuclear reactions because fingering convection does not mix the 
stellar matter deep enough. However, for an accretion of 1 M® 
(cyan dot-dashed line) beryllium is already reduced by a factor 5. 
If the beryllium abundance determination of the two components 
of the 16 Cygni system is confirmed, it may lead to a precise de¬ 
termination of the accreted mass. 



Fig. 10. Beryllium abundance profiles after the accretion of different 
masses: no accretion (blue solid line), 0.66 M e (red dotted line), and 
1 M e (cyan dot-dashed line) at the beginning of the MS. 


6. Summary and discussion 

The 16 Cygni system is particularly interesting for comparative 
studies of stars with and without planets. The two main stars of 
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this system have been observed in several ways, leading to very 
precise constraints. 

The most striking feature is the lithium abundance, which is 
slightly higher than solar for 16 Cygni A and lower by a facto r 
of at least 4.7 in 16 Cygni B dFriel etal.ll 19931 : Iking et al.lll997l) . 
The lithium value given for this later star is an upper limit, which 
means that it could be completely lithium depleted. 

We here presented detailed computations of these two stars, 
leading to models that precisely fit their asteroseismic frequency 
determinations, as well as their radii, luminosities, effective tem¬ 
peratures, and gravities. These models were computed by fully 
taking into account atomic diffusion of helium and heavy ele¬ 
ments. Their characteristics are given in Table Efl 

Then we discussed the importance of the lithium observa¬ 
tions for the two stars. We first showed that the very large lithium 
depletion observed in 16 Cygni B cannot be accounted for by 
the classical means of rotation-induced mixing or similar types 
of extra mixing. The observable parameter difference between 
them is too small to explain such a large difference in lithium 
as observed. We suggest that the accretion of heavy matter onto 
the planet host star 16 Cygni B in its early main-sequence phase 
induced a special kind of mixing, namely fingering convection, 
which led to a large lithium destruction on a relativ ely small 
timescale (cf. lVauclairll2004l : lTheado & Vauclairll20l 2l) . 

We recall that when a star accretes heavy matter onto lighter 
one, the induced //-gradient leads to a specific kind of hydro- 
dynamical instability, called fingering convection or thermoha¬ 
line convection, which mixes all the accreted matter downwards. 
Contrary to what is often assumed, the accreted heavy elements 
do not leave any signature in the stellar outer layers because of 
this mixing. On the other hand, the same mixing may transport 
lithium down to the layers where it is destroyed by nuclear reac¬ 
tions. 

We suggest that during the evolution of the planetary disk, 
at the beginning of the main-sequence phase, 16 Cygni B may 
have swallowed a fraction of an Earth-like planet that has led to 
fingering convection below the convective zone. This short-time 
efficient extra-mixing allowed the transport of lithium down to 
the nuclear destruction region. We have seen that two thirds of 
an Earth-mass planet would be enough to account for the ob¬ 
served upper value of lithium in 16 Cyg B, as shown in Fig. [7] 
Obser vatio ns of beryllium have been reported in these two stars 
bv IDelivannis et al.l ( 2000 ). They found no difference in beryl¬ 
lium between the two stars. If confirmed, this means that the 
accretion of planetary matter could not have exceeded one Earth 
mass. 

We also computed the detailed variations of the heavy ele¬ 
ment abundances due to atomic diffusion along the two evolu¬ 
tionary tracks leading to the best models for 16 Cygni A and B. 
Atomic diffusion alone leads to carbon depletion by about 0.05 
dex and to magnesium depletion by about 0.03 dex. The result¬ 
ing abundances in th e two stars are very similar, as observed by 
ISchuler et al.l (1201 ll) . Differences as suggested bv lRamfrez et akl 

1 Note added in proof: After this paper was accepted, Travis Met¬ 
calfe draw our attention to his new article accepted for publica¬ 
tion in ApJ "Asteroseismic modeling of 16 Cyg A & B using the 
complete Kepler data set" by T.S. Metcalfe, O.L. Creevey & G.R. 
Davies, arXiv : 1508. 00946v2. The authors revise the data given in 
iMetcalfe et all (120121') . particularly the age (7.0 ± 0.3 Gyr) and the com¬ 
position (Z = 0.021 ± 0.002, Yi = 0.25 ± 0.01). They acknowledge that 
the helium content they derive may be slightly underestimated due to 
their neglecting atomic diffusion of heavy elements in the models. It is 
interesting to stress that we indeed find a slightly higher helium content 
and a younger age than they. 


120111 and lTucci Maia et al.l20 14l cannot be explained in this con¬ 
text. 

In summary, the 16 Cygni system gives evidence of fingering 
convection induced by the accretion of planetary matter, which 
is probably the unknown process invoked in the past by several 
authors to account for the large lithium differences between the 
two stars. 

This result can be generalized to all planet-hosting stars, 
which may accrete planetary matter in a rand om wa y, w ith 
various accretion rates. As suggested by [Theado et al , 1 d2012t) . 
this could lead to an average lithium abundance that is lower 
in planet-host i ng sta rs t han in other stars, as cla imed by 
llsraelian et al.l d2009l) and iDelgado Mena et all d2014l) . In this 
framework, precise lithium abundances in seismically observed 
stars with or without observed planets would be very useful. Fur¬ 
thermore, beryllium observations in the same stars would help 
constraining the mass accreted onto the star. 
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